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Non-local f{R) gravity was proposed as a powerfull alternative to general relativity 
(GR) . This theory has potentially adverse implications for infrared (IR) regime as well 
as ultraviolent(UV) early epochs. However, there are a lot of powerful features, making 
it really user-friendly. A scalar-tensor frame comprising two auxiliary scalar fields, used 
to reduce complex action. However this is not the case for the modification complex 
which plays a distinct role in modified theories for gravity. In this work, we study the 
dynamics of a static, spherically symmetric object. The interior region of spacetime 
had rapidly filled the perfect fluid. However, it is possible to derive a physically based 
model which relates interior metric to non-local f{R). The Tolman-Oppenheimer-Volkoff 
(TOV) equations would be a set of first order differential equations from which we can 
deduce all mathematical (physical) truths and derive all dynamical objects. This set of 
dynamical equations govern pressure p, density p, mass m and auxiliary fields {-0, ^}. The 
full conditional solutions are evaluated and inverted numerically to obtain exact forms of 
the compact stars Her X-1, SAX J 1808.4-3658 and 4U 1820-30 for non-local Starobinsky 
model of /(□“^R) = . The program solves the differential equations 

numerically using adaptive Gaussian quadrature. An ascription of correctness is supposed 
to be an empirical equation of state ^ = a(l — e ^ pc) for star which is informative in 
so far as it excludes an alternative non local approach to compact star formation. This 
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model is most suited for astrophysical observation. 

Keywords: Higher-dimensional gravity and other theories of gravity; neutron stars; ther¬ 
modynamic processes; conduction; convection; equations of state. 

PACS numbers: 04.50.-h,97.60.Jd,95.30.Tg 


1. Introduction 

Modern cosmology is based on the relativistic description of large scale Universe. 
The situation is made worse when the Universe contains expansion and so acceler¬ 
ating. We are restricted to obey different observational data in favor of an acceler¬ 
ating universe at very large scaled Accelerating expansion is largely responsible 
for the fact that each galaxy in large scale are expanding the capacity to expand 
themselves. They provide evidence of the accelerating expansion of Universe and 
hence for the existence of the mysterious dark energy which drives this process. 
There is no doubt that these data are modifying the classical description of gravi¬ 
tation, general relativity (GR). We will engage other theories in modifying our GR 
to finally understand this epoch. Universe may be accelerating due to matter and 
geometry of the full action. Two contrasting methods for accelerating scenario have 
been identified. At the first candidate, an exotic matter fluid was used in the full 
action of theory to realize a variety of accelerating and decelerating phases. They 
had to prove that GR respond to cosmological needs in a world of accelerating ex¬ 
pansion that we observed earlier. At the next stage, a geometrical modification was 
used in the action to study a variety of accelerating and decelerating epoch^^US|). 
As individuals, there are things we can perform to improve the situation, firstly as 
geometrical modifications and secondly as fluids. 

Firstly there is the belief that the Buchdahl was the founder of the f{R) mod¬ 
ified gravity.E^ A model is proposed to explore the viable range of parameters of 
cosmological background due to replacing the Ricci scalar R by an arbitrary func¬ 
tion f{R). Proposed action by nonlinear higher terms should be applied to the large 
scale. Among the proposed models it has been considered that Universe expansion 
caused in the late time where low curvature corrections were formed as space time 
were filled. Similar difficulties may apply to the higher order corrections as Gauss- 
Bonnet (GB) term, G = R^ — . Indeed, in four dimensional 

Einstein-Hilbert action, these GB terms have no contribution if they appear in 
non-minimally coupled form. The reason backs to the topology and topological in¬ 
variants. Coupling of GB term to matter fields provides more dynamical features 
than before. In comparison to f{R)^ the modified GB gravity was proposed to be 
the /(G) gravity as model for dark energy in.l^ Several interesting cosmological 
features of this type modified gravity were studiecP^.^ One of the most impor¬ 
tant motivations of GB gravity is that it arises naturally from string theory. Also, 
it appears as a non-minimally coupling term to higher order scalar fields The 
list of possible modifications is not limited upto the above mentioned models. The 
non-local models are also so interesting theories and enough popular as modified 
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gravities. This last case, is studied in our paper. We’ll present a short review of the 
subject in Sec. 

These mass densities that are localized then immediately collapse, releasing ther¬ 
modynamics that can cause a hydrodynamic pressure impose on the mass distribu¬ 
tion. Following the collapse of the matter distribution, mass transformed itself into 
a compact star. The compact star has a finite size and a huge amount of mass.^^ It 
was an open competition between different massive objects that enabled the neutron 
star carriers the opportunity to exist in the space time. The neutron star has very 
strong surface gravity and thermodynamics. We can detect these massive objects 
by looking at the Doppler shift in spectral lines emitted by atoms in the surface. 
The major elements mass M ~ 1.4M0 and radius R ^ lOKm of inputs are used to 
control the whole dynamics. Specific surface gravity known as the relative surface 
gravity of the neutron star also may be studied. Not much gravity on the earth, 
you can find more on the compact star k = 2 x 10^^! The field of force associated 
with compact star have both electric and magnetic components and contains a def¬ 
inite amount of electromagnetic energy. This means there is a greater difference 
between individual plants which makes the species and neutron stars. For some 
basic parameters, for example emission above-ground level, electromagnetic effects 
were stronger than thermodynamic effects. 

It was shown that the thermodynamic parameters, pressure p, energy density 
p, the mass function M and radius R are related according to a set of ’’state” 
equations, called as TOV equations in GRP^.^^It defines a set of first order differ¬ 
ential equations for a spherically symmetric metric which is filled with the perfect 
fluid with pressure and energy density. One main goal of our paper is to derive 
explicitly, the forms of TOV equations for non-local f{R) gravity. Recently TOV 
equations and dynamics of stars have been investigated for different types of modi¬ 
fied gravity models from f{R), /(G) and f(T) (T is torsion ) numerically and in a 
non-perturbative schem^^.l^ 

Our plan in this letter is the following scheme: In Sec. © we present non-local 
f{R) gravity as an alternative theory for gravity. In Sec. ([3]) we derive equations of 
motion for a spherically symmetric star. In Sec. (|4]) we pass to the dimensionless 
parameters and we redefine all functions to obtain TOV equations. In Sec. m we 
study an isotropic model of compact star using astrophysical data. We conclude and 
summarize in Sec. ®. Some preliminary formulae are presented in Sec. (??). 


2. Non-local f{R) gravity 

Non-local corrections to the Einstein-Hilbert action proposed an attempt to obtain 
a ’’healthy” version of GR with added quantum loop corrections^.®! The simple 
problem was how to explain the current acceleration expansion of the Universe 
and to get the large numbers from inverse differential operator (s). Indeed, this 
modification is refereed as an IR non-local modification of General Relativity. 
After the original one, another model proposed as non-local F{R) gravitj®!-.®! It 
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was essentially a viable IR modification of the original f{R) gravity. Our study will 
be started from this motivated idea (for a dicussion on singularities and cosmological 
aspects se^^. Let us start by the appropriate form of action for non-local f{R) 
gravity: 

5 = y'd4:rv^|^R(l + /(n-iR))+£ 

matter } • ( 1 ) 

Here / is an arbitrary function of i? , □ = di}j stands for 

Laplace-Beltrami (dAlembertian) operator in equations of motion. We also adopted 

the commonly used signature of the metric as (-1-). With this signature, 

the curvelinear derivative operator and the Riemann tensor read as the following: 


- P^.Pa, (2) 

r)(y _ pi per _ o per . pu; per _ pcj per /o\ 

^ fJ.l'P ^ pp-L ojjy ^ pi/^ ojp \'^) 

The above action m may be recast to the following scalar-tensor form using a pair 
of auxiliary (may be unphysical) scalar fields 


S= d^xyZ—g 


2«;2 


{R{l + f{^))+^{ai;-R)} + C 

matter 


(4) 


= / d^Xyf^ 


^ {-R (1 + /(V')) - - iR} + ^matter 


The equations of motion are supplemented by a set of Euler-Lagrange equations for 
^ in the following form: 

^ = 0, U^ = R (5) 

The above equation may be recast to the following form ?/; = If we substitute 

this equation into (j4]), we reobtain (HD- Equations of motion for metric tensor g^^ 
is obtained from -4^ = 0: 

Sgij.v 

0 = {R (1 + /(V') - e) - - R^. (1 + /(V^) - 0 (6) 

+ l (5^C5,V’ + - V^V,) (/(V-) - 0 + ■ 

If we write the equation of motion for tp we get the following equation of motion: 


0 = + f\P^)R . 


(7) 


Our aim here is to derive explicit forms of dSlIzl) for spherically symmetric static 
conhguration of compact, neutron, quark E^tars. 


'Quark (q) is a fundamental fermion that has strong interactions 
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3. Spherically symmetric model of compact stars 

We suppose there must be a broadly believable compact star in static-spherically 
symmetric coordinates given by system = (ct, r, 0, (p) in the following form: 

ds^ = - e^^dr^ - r^{d9^ + sin^ 9d<p^). (8) 

We also assume that the matter fields are same which have been used in the interior 
parts and are in the comoving motion. The appropriate energy-momentum tensor is 
expressed as = diag{pc^, —p, —p, —p). To keep the homogeneity, we assume that 
^ = ^(r), tit = 'tp{r). We insert ([5]) in dH]) and using the formula given in Sec. (??), 
we obtain the diagonal components of ([HI) for (/i, = {ct,ct), and = (r, r) 

are given by the following differential equations: 

1 2A' 1 — / \ 

tt: --^V + (l + /(^)-e)(-—+ —(9) 

- - A') {tp'fp - ^') = 


rr : {I + /(</)) - C)(— + ^ 2 ) + ^ “ (“ + “ O = -K^e^^p{10) 

r 2 r 


For scalar field il) we rewrite (O: 

psi : V'" + f- + 7 - aOV'' - 7' - + 7A' - -(7 - A') 

\r / r 


1 — e 


2A 


-A' 1- 


o 2 A 


+ 


= 0 ( 11 ) 


) =112) 


and similarly for © we obtain: 

xi: (^+^'-A')e' + 2/v,(<)>" + 7"-7A' + : 

The trace of the equation of motion ([6]) provides another useful equation: 

-R (1 + fW - C - - 3n(/(V’) - C = -K'^ipc^ - P)- (13) 

We are licensed to write the equation on to metric (|5]) 

Trace : 2(</>" + - 0'A' + -(</>' - A') +-^) (1 + /(V-) - C - 7^' 

\ r / 


r 

-3 (■*/'"/</- + - 0 + (^ + 7 - A')(V'7'0 - ■?') = -p). (14) 


4. Tolman-Oppenheimer-Volkoff Equations 

The gravitational equations of motion must be supported by an appropriate hydro¬ 
static equation for the matter fields inside the star. This equation is nothing just 
the familiar continuity equation for energy-momentum tensor : 
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If we put V = r in this equation with the metric (|5]), we get: 

^ =-{P + ■ (15) 

dr 

Now, keeping v = t in mind, a question may arise that this completely vanishes the 
possibility of constructing new hydrodynamic equation. We will note that equations 
(I9ll()lllll2ll4ll are reduced to GR by any f{R) = R model in the action. 

We need to find the forms of the TOV equations. It is safe to replace the metric 
part A of the equations with a mass function M = M{r): 

2GM GdM 1 


,-2A 


= I - 


1 - e"2^(l - 2rA') 


(16) 


c^r c^dr 2 

Now, we must rewrite (I9I10I11I12I14I) in terms of {^, ^,p}. Also, it is adequate 
to write equations in the dimensionless forms. It is desired to introduce the next set 
of the dimensionless parameters for field equations, 

pMqC^ 


M 

GMp, 


mM( 


©) 


rgr, p 


pMq 

3 ’ ^ 


(17) 


where rg = = 1.47473ArM and Mq Stands out for its mass of the central 

Sun. Using these dimensionless parameters we rewrite the continuity equation as 
the following: 


P' 


(P + P)' 


But (fTBl) converts to the following: 


A' = 


1 - 


dm 
, dr 


^2 2m 
r 

Furthermore, we have metric function as an equation in terms of the mass: 
6^^ = 


(i-—r‘. 

1 - = (l - 

V r / 

V 2m) 


(18) 


(19) 


( 20 ) 


A very widely expressed system, TOV equations, connecting metric to matter is 
also found in the following forms: 

dm 


tt : 


+ (1 + /(</,) -^)(-1 mar 


1 


2(P + p) 


3 2m _ 2 ^ 
r 

ml — 




) 


r dm 
m dr 


p2 2m 


zm _ ^ 




: (1+ /(</>)-0(- 


2p' 


ip + p) 






r p + p 


)«.'/.-f')= 


^ _ zm 


2m ’ 




2{p + p) 


( 22 ) 
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psi 


i: r+(- 

\r 


p' 


1 r dm , , f . I . / . o 


p + p 


r p + p ^ — 17'^ ' \p + pJ Vp + py 

dm r\ 1 r dm -i 

= 0. (23) 


/I r dm n f r dm 

rn P , 2, P rn l--^ 


_ m dr I / 

^ _ 1 


p + p r 


.2 2m 


-1 


)- 


’^^(1-2^) 


Similarly, an equation could find in ([7]) : 

'2 p' 


i: e+(- 
v? 

■(—)' 

Vp + p/ 


'-.r p + p r 

/ V 9 /I r dm o , / 

_L ^ P _ 2/ P 

Vp + p 


- JL am / n' \ I 

( t 1 "P dm ^ 

P , rn 1 - 

n + 0 ^ — 1 y' 


p + p 2^ — 1 rVp + p r 


ifl- 

2m y 


Trace 


: 2f- f^)'+f^)\^- — 
V Vp + p/ Vp + p/ p-\-pr^ ^ 

1 


^ _ _?2_ 

m dr 


c\ f IP dm 

2/ P , , 

\ ..9 9m. -I / 


r'p + P 


- 1 

(i+/(V')-e)-V'r 


(24) 


= 0 . 


(25) 


-3 


{r u+v^u^-£.") +V- 


2m J 

2 p' ml - ^ 


27r^)(V''^-n 


V p + p ^ — 1 

Try adding up our model f{R), we solve (I21I22I23I24I25I1 numerically. 


87r(p - p) 


5. An empirical model for the astrophysical objects 

In previous section we derived the full set of TOV equations for a generic model 
of non-local f{R) gravity. Our aim in this section is to investigate a simple model 
of compact stars. By compact star, we mean some relativistic massive objects with 
tiny size and high density.!^ These astronomical objects have the mass of order 
M ~ Mq and the radius R ^ 10AM. For three types of the astronomical candidates 
the metric functions {A, tj)} were obtained as simple quadratic functions of radial 
coordinate r!^ 


2X = Ay, 2(t> = By+C, (26) 

where A, B and C are physical constants to be evaluated using astronomical data 
given in the table I. In GR, compact stars with these metric functions were studied. 
In modified gravity from f{R) to the GB model^^,^ there are some interesting 
features. In reference,l^a model for neutron star constructed numerically for a class 
of viable models of /(G) gravity. Our aim in this letter is to investigate physical 
properties of compact stars in this non local f(R) gravity (se^^ for some recently 
reported works). We will study the dynamical stability, energy conditions and red 
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shift properties of an isotropic model of compact star. The model which we address 
here is the non-local form of the one proposed by Starobinsky for inflation!^ 

f{a-^R) = D-^R + a(p-^Ry. (27) 

Even the best known neutron stars often have great uncertainties in their masses 
and radii. In fact no neutron star has a really precise radius measurement, and it 
is rare to have even rough measurements of both radius and mass for one star. 
The mass of strange star, describes a fit with two minima of error, one is O.SM©, 
the other being 1.8Mq which is far more believable as it is consistent with good 
measurements of other systems,and both mass values are highly dependent on the 
details of this very complicated model, in contrast with the simpler methods that 
produce good mass measurements in other systems. This binary is complicated by 
x-ray heating of the primary, among other things. Since even the worli^ say that 
the O.SMq measurement cannot be favored by their analysis, it is unreasonable 
to take it as definitive. The R ^ 7Km radius for SAX JI808 may extract fronP^ 
which it seems to be still popular in the particle literature although the astronomy 
literature shows it to be wrong, as iiP^ and others. We still cannot figure out how 
can we obtain the mass and radii for the star 4U 1820-30 (for a recent work see for 
exampl^^ for a detailed and careful discussion, recent analysis including statistical 
and (numerous) systematic sources of error). 


Strange star candidate 

A{km ^) 

B{km ^) 

Her X-1 

0.006906276428 

0.004267364618 

SAX J 1808.4-3658 

0.01823156974 

0.01488011569 

4U 1820-30 

0.01090644119 

0.009880952381 


• Exact solutions for auxiliary fields 

are completely integrable : 

Qrrdr 


'0, ^ : the equations of motion dSEl) 




\/- 9 {y)R{y)i.r - y)dy + C 2 , (28) 


here R{y) is (HD) . For ^ using ([7]) we obtain: 

'J-g{y){r - y)f^{y)n'ip[y)dy -b (74, (29) 


Where □^/'(y) was obtained from (14211 . However, it is their hard work, in¬ 
tegrating and simplifying to achieve their fields solutions (l28l [2^ . Fortu¬ 
nately, exact solutions (128112^ have given the graphically monotonic forms. 
In Fig. o, we observe that ^/’(r) never falls into a linear plot, which remains 
the polynomial. So auxiliary field '(/’(r) is to be the perfect polynomial plot 
against the radius. An example Her X-1 plot of field x(r) in an interior 
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r(Km) 


Fig. 1. Numerical plot of the 'tp{r) II28I I 


area is shown in Fig. ©. We conclude that ^(r) is always decreasing or re¬ 
maining constant and never increasing, so both of ^(r) are monotonic- 
increasing (decreasing) functions. For Starobinsky model, this behavior is 
a functional dependency of auxiliary fields on scalar curvature R of met¬ 
ric (ESI)- You may increase or decrease the Ricci scalar R m simply by 
extending or decreasing the {A,B,C}, or by increasing or decreasing the 
amount of fields tp, However, the increasing value of R points to the grow¬ 
ing appearance of second derivative of the scalar fields {□^,□'0}. At the 
same time both field(s) and growing Ricci scalar are increasing (decreasing) 
behavior. 

• Stability-conditions: with applying external radial perturbations at the 
inflow boundary, large radial structures develop naturally in the flow field 
due to sound effects. The velocity of sound drift, produced by spatially 
uniform perturbations is obtained by: 




dp 

dp 


We plot (1^ in Fig. ([3|). The probability (velocity) for the radial pertur¬ 
bations will peak at order ten to the one. Generally, however, there are 
instabilities with either pressure radial perturbations or with density to as- 
trophysical groups SAX J 1808.4-3658,Her X-1 and 4U 1820-30. One may 
increase or decrease the simply by extending or decreasing the pressure, 
or by increasing or decreasing the amount of density. For SAX J 1808.4- 
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Fig. 2. Numerical reconstructed plot of ^(r) using II29I I 

3658, numerical analysis showed a significant similar relationship between 
increasing proportions of sample SAX J 1808.4-3658, Her X-1 and 4U 1820- 
30. Initially, increasing the pressure perturbation to the density perturba¬ 
tion increases the velocity of sample. The plot may have a singularity, or 
the perturbation scheme may be inappropriate. But an instability ought to 
come after the singular point. The earlier mentioned sample SAX J 1808.4- 
3658 of astrophysical object is a such indicative case. However, looking back 
we should have instabilities in SAX J 1808.4-3658 a lot earlier than Her X-1 
and 4U 1820-30. Thermal instability of the stars with large perturbations 
is a function of the empirical parameters {A, H, C} of the sample. 

• Surface Redshift: we’ve an increase in the wavelength of radiation emit¬ 
ted by a celestial body as a consequence of the gravitational field. The 
gravitational redshift z of thermal spectrum detected at infinity can be 
computed as 

z = e"'^-l. (30) 

We plot (1^ in Fig. (j4]). It is observed that the largest redshift occurs for 
an emitter at the center of the star. 

• Energy conditions: by looking at the field equations in modified gravi¬ 
ties we are often able to arrange a mutually acceptable energy density peff 
and pressure Peff- These energy conditions have altered our understanding 
of the range of conditions under which energy transfering onto a region 
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Fig. 3. 


The probability (velocity) for the radial perturbations 


dPjPc 

dp! Pc 



Fig. 4. Redshift of photons iIMt emitted from the center of the star as a function of r for three 
different candidates SAX J 1808.4-3658,Her X-1 and 4U 1820-30. 


















April 7, 2015 0:52 WSPC/INSTRUCTION FILE TOV-NONLOCAL- 

[JMPA 


12 D. MOMENI et al. 

occursPl^ 


NEC 

^ Peff Peff ^ 0- 

(31) 

WEC ^ 

Peff > 0 and Peff + Peff > 0. 

(32) 

SEC ^ 

=4> Peff + 3peff > 0 and Peff + Peff > 0. 

(33) 

DEC 

Peff > 0 and Peff ± Peff > 0. 

(34) 


The essence of energy conditions is purely geometrical simple and, as Hawk¬ 
ing said, self-evidentP However, we also believe that energy conditions in 
their pure essence are effective. Fig. ([S]) showing the changed in order of 
decrease of pressure ^ El- The resulting graph ([S]) shows the rapid decrease 
in — in the three different candidates SAX J 1808.4-3658, Her X-1 and 4U 
1820-30 and the more gradual decrease through the 4U 1820-30 and SAX J 
1808.4-3658. The resulting graph ([6]) shows the rapid decrease in scaled pres¬ 
sure ^ in three different candidates and the more gradual decrease through 
the 4U 1820-30 and SAX J 1808.4-3658s. All the pressures get vanished on 
the radius of star and none of the pressures is remained at all. For Her X-1, 
the pressure is vanished near the r ~ 8.25 and for SAX J 1808.4-3658, at 
r ^ 6.5 and 4U 1820-30 at r ^ 11. The precise distances for the vanished 
pressures can be seen in the.^^^^^^ The levels of radius found numerically 
are comparable to those found by astrophysical data.^^^^^^ These two meth¬ 
ods clearly produce comparable data. Information verified and indexed by 
data from the jllElE]) would be easily cross-referenced in WEC and SEC. 
The satisfactions of the all energy conditions will be independently verified 
numerically. Then again we empirically verified that our non local model 
for compact star was correct. The role of non locality in compact stars 
was therefore verified in saturation of the all energy conditions. However, 
despite strong interest in their usage, a lack of fundamental test data and 
verified structural non local theory guidance is inhibiting uptake. 

• An empirical equation of state: 

Fig. m shows the EoS of a compact star calculated numerically. The 
results suggest that the appropriate form of EoS is given by the following: 

^ = a(l-e"''^). (35) 

Tiny density increases the pressure size and gave value to the hardenability 
of star. We claimed that increasing density p ” will thwart the aspirations 
of future neutron star ”. Quark stars are also required to attend a generic 
EoS which emphasizes linear relation p = Ap + B. However, if the density 
is decreased, you would be required to attend the quark star EoS at the 
non local f{R) model. Normally, someone contemplating these low densities 


^In all of these graphs pc = -y^,Pc 
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Fig. 5. Scaled energy density — of the star as a function of r for three different candidates SAX 
J 1808.4-3658,Her X-1 and 4U 1820-30. 



Fig. 6. Scaled pressure — of the star as a function of r for three different candidates SAX J 
1808.4-3658,Her X-1 and 4U 1820-30. 
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Fiff. 9. Numerical plot of the -^ 

® Pc Pc 


p <C Pc would attend quark star EoS. All parameters are shown in Table 
2. Some parameters offer a selection from a sample to which you give an 
option. For some objects parameters, for example 4U 1820-30, high pressure 
effects are stronger than pressure effects in other samples. Finally, the EoS 
parameters must also be consistent with a linearity in the quark EoS which 
leads to the observed samples. These parameters are tested in both density 
and pressure initial conditions. We mention here that the EoS (1551) can be 
addressed as a generalized exponential Virial EoS!^ 

^=exp(sS.,A„p---‘), (36) 

where the coefficients Km and the virial coefficients Bn are related, ks is the 
Boltzmann constant and T temperature. There are significant differences 
in the EoS of (1551) versus exponential Virial EoS (1551) . These differences 
have been categorized in terms of pressure of the background pb Pc- For 
instance, for anomalous differences all the normal gases with exponential 
Virial EoS have an pb of zero. These differences arise because one or more 
compact objects in our case has a pb value or indeed should not be omit. 
However, it is debatable whether these structural differences make this EoS 
any easier for the astrophysical purposes to investigate the thermodynamic. 
The differences between (1351) and exponential Virial EoS (l36l) are: (l35l) is 
made from the gravitational field, whereas exponential Virial EoS (l36)) is 
made from the high-density fluid. In conclusion, it is of utmost importance 
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p/p. 


p b ^ 

Fig. 10. Empirically computed quation of state -p- = a(l — e pc). Here cJa-,<yb standard 
errors for data fitting. 


to reiterate the differences between the ((35l) EoS and the exponential Virial 
EoS (1551) . 


Astrophysical strange star candidate 

a 

(7a 

b 

(7b 

Her X-1 

0.45526 

0.0198 

0.0152 

0.00124 

SAX J 1808.4-3658 

1.2242 

0.0563 

0.0191 

0.00162 

4U 1820-30 

1.1255 

0.0408 

0.0243 

0.00163 


6. Summary and conclusion 

Einstein gravity is a gauge theory of gravity. It should be modified to have more 
effective predictions and implications for recently observational data. One of the 
most accepted modifications of Einstein gravity is f{R) gravity and its extensions. 
It is assumed that we gain more information about gravity if we replace R by an 
arbitrary function f{R). Several cosmological aspects of this type of modified gravity 
have been investigated in literature. Specially the late and early time evolution. Non¬ 
local terms, induced by quantum effects can be considered as non-local higher order 
corrections to Einstein gravity. It is reasonable to consider both scenarios in a same 
context, as non-local f{R) gravity, a scenario which we studied in this letter. We 
derived the equations of motion for this non-local theory using a pair of auxiliary 
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scalar fields. As a motivated idea, we studied stellar structure using the modified 
forms of TOV equations. We obtained the set of equations of motion for a star in 
non-local form of f{R) gravity. At the same time the TOV equations were recast 
and some reconstruction took place in the system, when a non-local correction was 
inserted. It is asserted here that the dynamic can adequately describe, explain or 
understand such a naive relationship from the perspective of compact stars. The 
full conditional solutions are evaluated and inverted numerically to obtain exact 
forms of the compact stars Her X-1,SAX J 1808.4-3658 and 4U 1820-30 for model 

of - 1 - . The program solves the differential equations 

numerically using adaptive Gaussian quadrature. An ascription of correctness is 
supposed to be an empirical equations of state = a(l — for star which is 

informative in so far as it excludes an alternative non local approach to compact star 
formation. The differences between (1351) and exponential Virial EoS (l36t are: (1351) 
is made from the gravitational field, whereas exponential Virial EoS (I36p is made 
from the high-density fluid. This model is most suited for astrophysical observation. 
A theoretical perspective proposed by us, TOV equations for non-local f{R) theory, 
is helpful for understanding these non-local effects. 
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7. Appendices 

In this appendix we present different geometrical quantities which have been used 
in this letter. For metric ([ 8 ]) the following nonzero components of the symmetric 
connection are obtained: 


rL = </>', r22 = A', f 

r 44 = —r sin^ r 23 =-, FI 4 = — sin 0 cos 0 , r 24 =-, 

r r 

So, the nonzero (tt), (rr) components of the Ricci tensor read: 

R,, = + cf'X' - ^), 

r 

9 \' 

_ ^’X' - —. 


, = -re-2^ (37) 
= cot 6. (38) 


(39) 

(40) 


The Ricci scalar is as the following: 


R = -2e“2^(<b"-k(5'2 



A') 


1 — 

—^2 )■ 


(41) 


r 
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The operator □ is given by: 

□A(r) = -e-2^ (a" - A') A') . (42) 
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